Orbital selective Mott transition in multi-band systems: 
slave-spin representation and dynamical mean-field theory 
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We examine whether the Mott transition of a half-filled, two-orbital Hubbard model with unequal 
bandwidths occurs simultaneously for both bands or whether it is a two-stage process in which the 
orbital with narrower bandwith localizes first (giving rise to an intermediate 'orbital-selective' Mott 
phase) . This question is addressed using both dynamical mean-field theory, and a representation of 
fermion operators in terms of slave quantum spins, followed by a mean-field approximation (similar 
in spirit to a Gutzwiller approximation). In the latter approach, the Mott transition is found to 
be orbital-selective for all values of the Coulomb exchange (Hund) coupling J when the bandwidth 
ratio is small, and only beyond a critical value of J when the bandwidth ratio is larger. Dynamical 
mean- field theory partially confirms these findings, but the intermediate phase at J = is found 
to differ from a conventional Mott insulator, with spectral weight extending down to arbitrary low 
energy. Finally, the orbital-selective Mott phase is found, at zero-temperature, to be unstable with 
respect to an inter-orbital hybridization V , and replaced at small V by a state with a large effective 
mass (and a low quasiparticle coherence scale) for the narrower band. 

PACS numbers: 71.30+h,71.10.Fd,71.27.+a 
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I. INTRODUCTION 

The Mott metal-insulator transition plays a central 
role in the physics of all strongly correlated electron ma- 
terials. At a qualitative level, localization of the elec- 
trons can occur when the kinetic energy gain (typically 
given by the bare bandwidth) is smaller than the cost in 
on-site repulsive Coulomb energy (U). In recent years, 
dynamical mean-field theory (DMFT)ii^ has provided 
a consistent theoretical framework which has advanced 
our understanding of this phenomenon^, in particular 
through the study of simplified models such as the one- 
band Hubbard model. 

In real materials however, such as transition metal ox- 
ides, several orbital components are involved. Crystal- 
field effects and the Coulomb exchange energy (J) affect 
the energy of on-site atomic states, which no longer de- 
pend only on the total local charge as in the orbitally 
degenerate case. Furthermore, the inter-site hopping am- 
plitudes can be significantly different for different orbital 
components (due e.g to their relative orientations). It 
is therefore essential to understand how these effects can 
affect the Mott transition, and whether qualitatively new 
effects are possible when the orbital degeneracy is lifted. 

Recently, this question has attracted a lot of at- 
tention. In their study of Ca2- a; Sr :! ;Ru04, Anisimov 
et alA suggested that a partial localization could take 
place, in which some orbital components (with broader 
bandwidth) are conducting, while others (with narrower 
bandwidth) are localized (see also Refpi). Following 
this proposal, several studies have been performed in 
the model context, with controversial resultsZ&SiiSiii. 
Liebscb2ii£ initially challenged the existence of such an 
"orbital-selective Mott transition" (OSMT), on the basis 
of DMFT calculations. Koga and coworker a 9 ! 11 , on the 



other hand, did find an OSMT within their DMFT calcu- 
lations, and suggested that a unique transition is found 
only if J = 0. A symmetry argument was put forward to 
explain this finding. 

In this paper, a clarification of this problem is at- 
tempted, using DMFT and another, complementary ap- 
proach. The latter is based on a representation of fermion 
operators in terms of slave quantum spins, specifically 
forged to address multi-orbital models (Sec. IHIfl . A 
mean-field approximation based on this representation, 
similar in spirit to the Gutzwiller approximation, pro- 
vides a fast and efficient method in order to investi- 
gate the Mott transition in a wide range of parameters 
(Sec. llVjl . In Sec[Vl a detailed study of the previously 
unexplored regime in which one of the bands is much nar- 
rower than the other and J = is presented, using exact 
diagonalizations and quantum Monte Carlo methods in 
the DMFT framework. Finally (SecEB, the effect of an 
inter-band hybridisation is considered. 



II. MODEL 

The model considered in this paper is a tight-binding 
model for two bands, coupled by local interactions. The 
Hamiltonian reads H = Hq + Hi nt , where Hq is the non- 
interacting part: 
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in which d\ ma [di ma ) creates (annihilates) an electron on 
the site i, in the orbital m, with spin a. The e TO 's are 
crystal-field levels and fj, is the chemical potential, kept 
here for generality. In most of the paper however, we 
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shall focus on the case of zero crystal-field splitting (ei = 
£2 = 0) and half-filling of each band (i.e one electron per 
site in each orbital, which corresponds to fi = given our 
normalization of the interaction term). At the end of the 
paper, we shall also consider the possibility of a non-zero 
inter-orbital hybridization. 

The full interaction, in the case of degenerate bands in 
a cubic environmen1ji2*i£ reads: 
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where hi„ lcr = rii ma — 1/2. Following Castellani et ali£, 
the reduction of the Coulomb interorbital Coulomb inter- 
actions U' as compared to the interorbital U is related 
to the Hund's coupling J by: 



U' = U-2J 



(3) 



In the case of vanishing Hund's rule coupling J = the 
interaction vertex (= U(ni + n 2 — 2) 2 /2) thus depends 
only on the total charge, while if J ^ the inter-orbital 
interaction is weaker than the intraorbital one and be- 
comes sensitive to the spin configuration. 



III. SLAVE-SPIN MEAN-FIELD THEORY 
A. Slave-spin representation 

In this section, we introduce a new representation of 
fcrmion operators in terms of constrained ( "slave" ) aux- 
iliary fields, which proves to be particularly convenient in 
order to study the multi-orbital hamiltonian above. The 
main idea at the root of any slave- variable representation 
is to enlarge the Hilbert space and to impose a local con- 
straint in order to eliminate the unphysical states. When 
the constraint is treated on average, a mean-field approx- 
imation is obtained. Different slave-field representations 
will lead to different mean-field theories. The quality of 
the mean-field approximation can be improved by tai- 
loring the choice of slave fields to the specific problem 
under consideration. In general, a compromise has to be 
found between the simplicity of the representation, the 
number of unphysical states which are introduced and 
the possibility of an analytical treatment of the resulting 
me an- field theory. 

For finite-U Hubbard models, Kotliar and 
Ruckensteini^ have introduced a slave-boson repre- 
sentation which can be used in the present context, 
when appropriately generalised to multi-orbital mod- 
els (in the spirit of the Gutzwiller approximation 15 ). 



However, this method introduces many variational 
parameters. On the opposite, S.Florens and one of the 
authors introduced a very economical representation of 
the N-orbital Hubbard model with SU(N) symmetry 
based on a single slave variable, taken to be the phase 
conjugate to the total charge on a given lattice site 
(slave rotor representation)A§ii£ . However, this method 
is not appropriate when the orbital symmetry is broken, 
as in the present work. 

Here, we introduce a new slave-variable representa- 
tion^ especially suited for dealing with multiband mod- 
els, and addressing orbital-dependent properties. The 
basic observation behind this scheme is that the two pos- 
sible occupancies of a spinless fcrmion on a given site, 
rid = and 7i<j = 1, can be viewed as the two possible 
states of a spin-1/2 variable, S z = —1/2 and S z = +1/2. 
This representation has been widely used in the case of 
hard-core bosons. In the fermionic context however, one 
needs to insure anticommutation properties, and this is 
done by introducing an auxiliary fermion /, with the ad- 
ditional local constraint: 



Pf = s z 



(4) 



In this manner, one obtains a faithful representation of 
the Hilbert space, which reads: 
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|0> = | n/ = 0,5« = -l/2) (5) 
= d^\0) = \n f = l,S z = +l/2) (6) 



This constraint eliminates the two unphysical states 
\n f = 0, S z = +1/2) and \n f = 1, S z = -1/2). This rep- 
resentation is easily extended to the multi-orbital case, 
by treating each orbital and spin species in this manner. 
Hence a set of 2N spin-1/2 variables and auxiliary 
fermions f ma are introduced (m = 1, • • • , N is the num- 
ber of orbitals), obeying the local constraint on each site: 
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This constraint can e.g be imposed with Lagrange mul- 
tipliers fields Xima- (t). 

We now explain how to rewrite the original hamil- 
tonian illll'l) in terms of the slave spins and auxiliary 
fermions. We consider first for simplicity the case J = 0, 
since the case J ^ requires an additional approxima- 
tion, as discussed later. For J = 0, the interaction in- 
volves only the total electron charge on a given site, and 
therefore reads: 



H 



J=0 — 




In order to express the non-interacting part of the hamil- 
tonian, we need to choose an appropriate representation 
of the creation operator of a physical electron, d\ . 
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There is some freedom associated with this, since dif- 
ferent operators in the enlarged Hilbert space spanned 
by the slave-spin and auxiliary fermions can have the 
same action on the physical (constrained) Hilbert space. 
We have not used the obvious possibility d> — ► S + f\ 
d — ► S~ f. This representation is correct in the physical 
Hilbert space (i.e when the constraint is treated exactly), 
but it can be shown that additional mean-field approxi- 
mations based on this representation will ultimately lead 
to a problem with spectral weight conservation because 
S + and S~ do not commute. Instead, we have chosen the 
representation $ — > 2S X f*, d — > 2S x f, which is identical 
to the previous one on the physical Hilbert space, and in- 
volves commuting slave-spin operators. With this choice, 
the non- interacting part of the hamiltonian reads: 



Ha — y^m^~] ^S x na .Sj ma .(fl na fj m<J + h.c) 

m <ij>,a 
( f m ~~ M) fima fima 



At this stage, no approximation has been made, provided 
the constraint is treated exactly. 
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renormalisation of the hopping: jt m — 4(5? 
The quasiparticle weight is associated with a different 
quantity, namely: 



Z n 



(13) 



Note that it depends in general on the orbital, a key fea- 
ture for the physics that we want to address with this 
technique. Both the renormalisation of the mass and 
the quasiparticle weight are self-consistently determined 
from the solution of the quantum-spin hamiltonian 11(J|I . 
which describes the charge dynamics. As clear from 113L 
metallic behaviour for orbital m corresponds to long- 
range order in Sj^, while Mott insulating behaviour cor- 
responds to (S m ) = 0. 

At this stage, the slave-spin degrees of freedom are still 
described by a quantum spin hamiltonian on the lattice, 
and we therefore make the additional approximation (iii) 
of treating this model on the level of a single-site mean- 
field. We thus have to solve-the single-site spin hamilto- 
nian: 



H s — 2h r , 
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7;) + H mt [S m<7 } (14) 



B. Mean- field approximation 

Approximations will now be introduced, which consists 
in three main steps: i) treating the constraint on average, 
using a static and site-independent Lagrange multiplier 
\ m<J ii) decoupling the auxiliary fermions and slave-spin 
degrees of freedom, and finally iii) treating the slave-spin 
hamiltonian in a single-site mean-field approach. This 
last step is quite independent of the two previous ones, 
and can be rather easily improved on, as done inil. 

After the first two steps, one obtains two effective 
hamiltonians: 



m <ij>,cr 
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(10) 



with H int [{S lm a}] = [//2Ei(Ew SU 2 for J — 0. In 
these expressions, tffl and J%[' are effective hoppings 
and slave-spin exchange constants which are determined 
from the following self-consistency equations: 

= ^m{Si ma S ' j ma ) (11) 
Jm^ = tmifimafjma + fjma fima) (12) 

The free fermion hamiltonian @ describes the quasipar- 
ticle degrees of freedom. Their effective mass is set by the 



in which the mean-field h Ti 
from: 



is determined self-consistently 
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where z is the coordination number of the lattice, 
equation can be combined with (jT2"|) to yield: 
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(15) 



This 
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In this expression, the fermionic expectation value is 
to be calculated with the quasiparticle hamiltonian 
Within this single-site mean-field approximation how- 
ever, the renormalisation of the hopping becomes identi- 
cal to the quasiparticle residue since {Sf ma S x ma ) factor- 
izes into (Sf ma ) 2 . As a result, the quasiparticle hamilto- 



nian reads: 



H 



eff 



E 



(z 



(17) 



with e km = —t m /z J2j,n.n(i) e_k ' tnc F° ur i er trans- 
form of the hopping. Equations (|13ll4lltjll7|) and the 

constraint equation (J7J) self-consistently determine the 
variational parameters h m , X m and Z m = A(Sf na ) 2 . They 
are the basic mean-field equations based on the slave- 
spin representation, which will be used below. Solving 
these equations requires to diagonalizc the single-site spin 
hamiltonian 114(1. corresponding to a A N x A N matrix. 

Let us finally discuss the case of a non-zero Hund's 
coupling J ^ 0. The first three terms in are easy 
to treat in the slave-spin formalism, since they involve 
only density-density interactions and are thus directly 
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expressed in terms of the Ising components of the slave 
spins. They read: 



constraint equation J7J) reads in this case: 



U' 



EE s ^f + ^E(E s L„f - £ E(E s 



(18) 

In contrast, the "spin-flip" and (intra-site) "pair- 
hopping" terms (last two terms in J5J) are more diffi- 
cult to deal with, since they involve both slave-spin and 
auxiliary fermions operators. As a result, four-fermion 
terms are introduced which require additional mean-field 
decouplings. For simplicity, we choose to mimic the ef- 
fect of these terms by replacing them by operators which 
have exactly the same effect on the slave-spin quantum 
numbers of the Hilbert space, namely: 

-JE [^n-[ S ill^il S i21 + 5 'ill 5 '- i lT' S 'i2T' S 'i21_ 

i 

~ J E \_ S ti^ S ti[ S a] S i2i + S ^ S £>i S ii-\ S iii ' (19) 

i 

Despite the fact that these terms connect the physical 
and unphysical parts of the Hilbert space (and therefore 
would strictly vanish if the constraint was implemented 
exactly), it is reasonable to expect that they will qual- 
itatively describe the physics of the spin-flip and pair- 
hopping terms when the constraint is treated on average, 
because their action on the slave-spin quantum numbers 
is the correct one. Hence, we shall use the following rep- 
resentation of the interacting part of the hamiltonian for 
,7^0: 

i rn,a 

+^5Z(JZ Simo) ~ "77 5Z(^Z Simo) 2 
i,m cr i,cr m 

~^E [^m^m^i^W + ^ai^nt^m^i2i 

i 



C. Benchmarks 

In this section, we perform some benchmarks of the 
slave-spin representation and mean-field theory 



1. Atomic limit (J=0) 

In the J = case we check explicitly that the atomic 
limit (i.e. t m = 0) of the degenerate N-band model with 
SU(2N) symmetry is correctly reproduced. Indeed our 
equations simplify drastically in this limit (t m = h m = 
0), leaving only the X rna — A to be determined. The 



uf{h — A) = Z 



2N 

- 1 E M Q Qe 

Q=0 



-0 m-(Q-N) 2 + \Q] 



(21) 



where np{e) is the Fermi function, Z = 

and Q is the total number of particles. Solving numeri- 
cally this equation for A leads to the correct "Coulomb 
staircase", as shown in Fig^ as long as T <§; U. At 
high temperatures the fact that we have imposed the 
constraints only in average limits the accuracy, but in 
practice T U is not a severe limitation. 



4 
3.5 

3 
2.5 

2 
1.5 

1 

0.5 








FIG. 1: (color online). Filling vs chemical potential for a two- 
orbitals impurity (atomic limit of a particle-hole symmetric 
Hubbard model), U = 2, /3 = 50: within the slave spin mean 
field (full line) and exact result (dashed line). The Coulomb 
staircase is correctly reproduced up to temperatures of order 
~ U 



2. N -orbital Hubbard model with SU(2N) symmetry and 
large-N limit 

Here, we apply the slave-spin mean-field approxima- 
tion to the N-orbital model (m — 1, • • • , N), in the case 
where all bands have the same hopping, with J = 0. The 
results for the quasi-particle weight as a function of U, 
at half-filling, are displayed on Fig. [21 A transition into 
a Mott phase is found for U > U C (N). The exact large 
N behaviour of U C (N) in the limit of infinite coordina- 
tion (DMFT) is knowni& to be linear in N, the slope 
being U c /N = 8|e|, where e = /"^ deD(e)e, This slope 
is correctly reproduced by the slave-spin mean-field ap- 
proximation, indicating that this approximation becomes 
more accurate as N is increased. 

One can actually calculate analytically the critical 
value of the coupling within this approach, for arbitrary 
N, by performing a perturbative expansion around the 
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FIG. 2: (color online). Quasiparticle weight, obtained from 
slave-spin mean-field for the N-orbital Hubbard model at half- 
filling (with, from left to right: N = 1,2,3,4). The non- 
interacting density of states is a semi-circle with half band- 
width D Inset: Dependence of the critical U on iV. The exact 
large- N behaviour is obtained. 

atomic limit for small h m . This yields: 

U c = 8(N+l)\e\ (22) 

which coincides with the numerical determination in the 
inset of Fig. |2] We also note that (|2*2l is precisely the 
result of the Gutzwiller (slave-boson) approximation in 
the multi-orbital case. 



D. Comparison with slave bosons and slave rotors 

As suggested by the fact they yield identical values of 
U Cl the slave-spin mean- field theory has many similar- 
ities with the Gutzwiller approximation (GA). In fact, 
as shown on Fig. [21 the whole dependence of Z on U is 
identical to that of the GA. 

The slave spin representation has several advantages 
over the slave-boson representations that can be used to 
formulate the GA. One advantage is that the number of 
variables is smaller: 2N spin- 1/2 degrees of freedoms in- 
stead of 2 N slave bosons (one associated to each state in 
the Hilbert space, in the absence of symmetries). An- 
other advantage is that the number of unphysical states 
is smaller than in slave boson representations, because 
the Hilbert space spanned by the 27V quantum spins is 
finite by construction, while the Hilbert space associated 
with the slave bosons in the absence of the constraint is 
an infinite-dimensional one. This might be useful in con- 
sidering finite-temperature properties and the entropy of 
the model. 

A similar remark applies when comparing the present 
slave-spin representation to the slave-rotor representa- 
tion recently developed by S. Florens and one of the 
authorsiSiii . This representation is specifically tailored 
to SU(2N) symmetric models, and is very economical 



since it introduces only one slave variable. Specifically, a 
(slave) quantum rotor and auxiliary fermions are intro- 
duced on each site such that: 

^ima ~ fima e ' ' ^ima = fima^ ' (23) 

where the phase is conjugate to the local charge, corre- 
sponding to the local constraint: 

Y.ULjim a -\) = Li (24) 

ma 

in which Li = 1/id/dQi is the conjugate momentum to 
the phase. It is clear from these expressions that there is 
a close similarity between the slave-spin and slave-rotor 
formalisms. Two important differences must be noted: 
(i) a single slave-rotor variable is introduced for all or- 
bitals and (ii) the Hilbert space of the unconstrained ro- 
tor is infinite-dimensional, containing an infinite tower 
of charge states \l) which are physical only for \l\ < N. 
As a result, mean-field approximations in which the con- 
straint is treated only on average are less accurate when 
the contribution of these unphysical charge states become 
sizeable. This is particularly true in the weak-coupling 
limit. On Fig. [21 we compare the slave-spin and slave- 
rotor result for the quasiparticle residue in the one-band 
case, in order to illustrate this effect. 

On the whole, slave-rotors and slave-spins offer two 
useful representations, the former being very economi- 
cal and well suited to situations in which only the total 
local charge is involved (e.g in Coulomb blockade prob- 
lems^), while the latter is well suited to the investigation 
of orbital-dependent properties, as in the present article. 
Both methods are easy to implement at a very low nu- 
merical cost, hence allowing for a fast and efficient inves- 
tigation of the phase diagram and phase transitions in a 
wide range of parameters. 

IV. ORBITAL-SELECTIVE MOTT 
TRANSITION WITH SLAVE-SPIN MEAN FIELD 
THEORY 

In this section, we use slave-spin mean-field theory in 
order to study the two-band model with unequal hop- 
pings. The non-interacting density of states of each band 
is taken to be a semi-circle (of half-width D\ = 2t\ and 
Z?2 = 2<2 < Di), corresponding to a Bethe lattice with 
infinite connectivity z — oo and nearest-neighbour hop- 
pings ti^/s/z. No crystal-field splitting is introduced 
(ei = £2 = 0) and we restrict ourselves to the case 
in which both bands are half-filled ((n\) = (712) = 1). 
The model is thus particle-hole symmetric, implying that 
the chemical potential \x = and Lagrange multipliers 
Ai = A2 = 0. The parameter space was explored for 
U > 0, J = -f- 0.5U i.e. U' = U - 2J = -f- U, and 
the ratio between the two bandwidths £2/^1 = -7- 1.0. 
For the study of the Mott transitions in this model we 
monitor the quasiparticle weights Z m — A(S^ na ) 2 . 
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FIG. 3: (color online). Quasiparticle weight of the one-band 
Hubbard model, obtained with slave rotors (thin line), slave 
spins, and the Gutzwiller (slave boson) approximations (thick 
lines). The latter two actually coincide. The small-17 be- 
haviour of the slave-rotor approache is due to the larger num- 
ber of unphysical states (see text). 



Fig. 0] displays the phase diagram within slave-spin 
mean-field theory for the bandwidth ratio t^/ti = 0.5. 
Three different phases are found: at small U both bands 
are metallic (i.e. Z m ^ 0), at large U both are insulat- 
ing (Z m = 0), and in between an "orbital selective Mott 
phase" (OSMP) is found in which only the band with 
largest bandwidth has Z\ ^ 0, while the narrower band 
has Z 2 = 0. 
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FIG. 4: (color online). Phase diagram (U vs U') for 
<2/ti = 0.5 at T = within the Slave Spins mean field theory. 
Inset, same diagram obtained with Exact diagonalization- 
Dynamical Mean Field Theory (ED-DMFT) irA The dotted 
line indicates J = 0, i.e. U — U' . In "phase I" both bands 
are metallic, in "phase II" both bands are insulating. "Phase 
III" is the orbital-selective Mott phase. 

In the inset of Fig. 01 we reproduce for comparison the 
result of Koga et al. 9 obtained within DMFT. Qualita- 
tively, one sees that the slave-spin mean-field compares 
rather well to the DMFT results. There are quantitative 
differences in the critical values of the couplings U and 



U' , a well known feature of Gutzwiller- like approxima- 
tions. Also, the linear dependence on U' of the upper 
boundary of the OSMP phase is due to the simplified 
treatment of the spin-flip and pair-hopping terms dis- 
cussed above (as indeed confirmed by the results of sec- 
tion EI3i. 

There is however one significant qualitative difference 
between the slave-spin results and those of Koga et al. 
(inset). We find that the endpoint of the OSMP phase 
does not lie exactly on the U — U' line. Hence, within 
the slave-spin mean field, the Mott transition becomes 
orbital-selective (OSMT) only when J exceeds a criti- 
cal value. This is a rather significant finding, since for 
,/ = the interacting part of the hamiltonian (Hi nt ) 
has full SU(A) spin-orbital symmetry, while for J ^ 
the symmetry is lower. In Ref£, it was argued that 
indeed the enhanced symmetry of the J = case pre- 
vents an orbital-selective Mott transition to occur. Our 
finding that a critical value of J is needed to induce an 
OSMT (for t 2 /ti = 0.5) suggests that symmetry consid- 
erations on Hi n t may not be essential to the existence 
of an orbital-selective transition. After all, the difference 
in bandwidths breaks the SU (4) symmetry from the ki- 
netic energy part of the hamiltonian. In order to study 
this issue in more detail, we perform in the next section 
a detailed study of the J = case. 



A. OSMT at J = 

In this section, we focus on the J — case, for which 
Hi nt has full SU (4) symmetry, and explore the nature of 
the Mott transition in the full range of bandwidth ratio 
from t^jti = to £2/^1 = 1. 

We find that the two bands undergo a common Mott 
transition at a single value of U — U c as long as the band- 
width ratio exceeds a critical threshold: t^/t\ > 0.2. In 
contrast, for £2/^1 < 0-2, an orbital- selective Mott phase 
is found, despite the enhanced symmetry of the interac- 
tion term. Fig. 0] displays our result for the phase dia- 
gram as a function of t%jt\ and U jt\. All transitions are 
found to be second-order when J — 0. On Fig. [5] the 
quasiparticle weights of each band is plotted as a func- 
tion of U, for several values of t-xjt\. The localization of 
the narrower band manifests itself as a kink in the quasi- 
particle weight of the wider band. As U is increased fur- 
ther, the wide band in turn undergoes a Mott transition. 
We observe that, within slave-spin mean-field, the quasi- 
particle weight of the wider band in the orbital-selective 
Mott phase coincides with that of a single-band model. 
This is because the slave-spin mean-field neglects charge 
fluctuations of the localised orbital, so that the physi- 
cal behavior of the wide band becomes effectively that of 
a one-band model as soon as the narrow band becomes 
localized. 

Our finding of an orbital-selective Mott transition at 
J = when £2/^1 is small enough, within slave-spin 
mean-field theory, raises two questions. First, is this find- 
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FIG. 5: Dependence of the critical U on the ratio tz/ti at 
J = and T = 0. All transitions are second-order. 



ing an artefact of the slave-spin approximation or does 
it survive a full DMFT treatment (i.e is it a genuine fea- 
ture of the infinite-coordination model ?). Second, docs 
this invalidate the argument based on the symmetry of 
H int ? The first question will be addressed in detail in 
Sec.0 in which a DMFT study will be performed, using 
exact diagonalization and Quantum Monte-Carlo tech- 
niques. We will show that indeed, a transition does exist 
at J = when ti/tx is small enough, but that the na- 
ture of the intermediate phase (OSMP) at low-energy 
is a rather subtle issue. In order to address the second 
question, let us briefly recall the symmetry argument of 
Koga et al&. The argument relies on the gap to charge 
excitations in the insulating phase. When J = 0, charge 
excitations mix the two orbitals because of the enhanced 
SU(4) symmetry. Instead, for J ^ 0, the charge exci- 
tations of lowest energy are independent in each orbital 
sector. As a result, it is reasonable to expect (at least 
when the kinetic energy term is treated in a perturbative 
manner) that the system can sustain two different charge 
gaps when J ^ while the gaps might coincide for J = 0. 
We observe however that this argument applies to the in- 
stability of the large-U Mott phase (in which both bands 
are gapped) when U is reduced, and suggests that, for 
J = 0, the Mott gap closes at the same value of U for 
both bands when U/t\ is reduced. It does not preclude 
however that a transition into an intermediate phase does 
exist, in which the "localised" band (with the narrower 
bandwith) is not fully gapped. As we shall find below, 
there is indeed clear evidence from the DMFT results 
that the orbital-selective Mott phase at J = is not a 
conventional Mott insulator and that the narrow ("lo- 
calised") band does have spectral weight down to zero- 
energy in this phase. Obviously the slave-spin mean-field 
approach is too rudimentary to be able to capture these 
fine low-energy aspects, but it is remarkable that it does 
allow us to infer correctly that an intermediate phase is 
indeed present. 
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FIG. 6: Z m at J = and T = for t 2 /ti = 0.5 (top), 0.25 
(middle), 0.15 (bottom). 



B. Dependence of OSMT on J 

Having clarified the situation for J = 0, we come back 
to the effect of a non-zero J, still within the slave-spin 
mean-field approximation. Fig. [7] shows how the phase 
diagram as a function of the bandwidth ratio ta/ii and 
of U/ti is modified for J ^ 0. One sees that a finite 
J enlarges the orbital-selective Mott phase and favors 
an OSMT. The critical ratio (ia/tijc below which an 
OSMP exists increases significantly, e.g. t%lt\ ~ 0.55 for 
J = 0.Q1U. With increasing J, (t2/ti) c tends towards 1. 
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Hence a common Mott transition for both bands is re- 
covered for all values of J only when t\ — t 2 . For a given 
bandwidth ratio t^/ti, the Mott transition is orbital- 
selective for J/U > (J/U) c . The dependence of this crit- 
ical ratio upon t%jt\ (he the location of the enpoint of 
the OSMP phase) is displayed in Fig. |H1 (It should be 
noted however that this critical ratio (J/U) c is underes- 
timated by our simplified treatment of the pair-hopping 
and spin-flip term). Finally, we found that for finite J, 
the insulator to OSMP transition remains second-order, 
while the metal to OSMP transition becomes first-order. 
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FIG. 7: (color online). Widening of the OSMT zone with 
increasing J/U at T = 0. Transition lines are shown for (from 
left to right) J/U = (dashed),0.01, 0.1 
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nian in which the spin-flip and the (on-site) inter-orbital 
pair-hopping terms are dropped, namely: 

H m t — U n im ^h lmi + (U - 2 J) 2J n ila h l2s 

im ia 

+ (U-3J)^n ila n i2a . (25) 

ia 

which is easily represented in terms of slave-spins as: 

H int = — y^y^ si na ) 2 + 

i m,cr 

+^E(E s ^) 2 -iE(E^) 2 ( 26 ) 

i,m a i,a m 

This study is also motivated by a comparison to Quan- 
tum Monte-Carlo treatments in which the spin-flip and 
pair-hopping terms are not easily treated. On Fig. El we 
display the slave-spin phase diagram found for t%/ti = 
0.5 at zero-temperature. On sees that the OSMP shrinks 
dramatically (albeit the two transitions do not actually 
merge). This finding sheds light on the results by Liebsch 
in refi 7 ' 8 ' 10 . Because this study was based on Quantum- 
Monte Carlo, hence neglecting the spin-flip and pair- 
hopping terms, it is natural that the orbital-selective 
phase can be found only in a very narrow range of pa- 
rameter space. The key role of spin-flip and inter-orbital 
pair-hopping terms for the OSMT was actually pointed 
out in the recent work of Koga et aliii. 

In Fig. we display the region in ti/t\, J/U parame- 
ter space where the Mott transition is found to be orbital- 
selective, analogously to Fig|SJ The critical ratio (J/U) c 
is found to be roughly exponential in t 2 /t\. In contrast 
to the case of the full hamiltonian (Fig[SJ), we find that no 
OSMT exists when ti/t\ exceeds a critical bandwidth ra- 
tio t 2 /ti ~ 0.6, for any J/U. Beyond this value no OSMT 
is found at any J/U. (Note that the upper critical line 
at large J/U corresponds however to the unphysical case 
of an attractive Coulomb interaction due to U — 3 J < 0). 

Finally, we emphasize that the orbital-selective Mott 
phase, which exists only in a very narrow range of cou- 
plings for the simplified interaction (|25|l at T = 0, is actu- 
ally enlarged at finite-temperature, as shown in Fig. 1111 



FIG. 8: (color online). Dependence of the critical J/U above 
which the Mott transition becomes orbital-selective, on the 
ratio t 2 /ti at T = 0. 



V. DYNAMICAL- MEAN FIELD THEORY FOR 
J = AND NATURE OF THE OSMT PHASE 



C. Role of the spin- flip and pair- hopping terms in 
the Hund Hamiltonian 

In order to clarify the role played by the different terms 
of the interaction ©, we have also studied the Hamilto- 



In this section, we study the two-band model with 
J = using dynamical mean-field theory. Our goal 
is to determine whether the transition into an orbital- 
selective Mott phase (OSMP) found within the slave- 
spin approximation at small enough t%/ti is indeed a ro- 
bust feature, and to shed light on the possible low-energy 
physics of this phase. Within DMFT, the lattice model is 
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FIG. 9: (color online). Phase diagram (U vs U') for ti/tj = 
0.5 at T = for a two-band Hubbard model without spin flip 
and pair hopping terms in the interaction. 
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FIG. 10: (color online). Critical J/U for the model without 
spin-flip and pair-hopping as a function of t-z /t\ . At small ra- 
tios of the bandwidths, the system displays an OSMT above 
a critical J/U ratio which vanishes at i2/ii = 0.2. In less 
anisotropic systems large Hund's coupling are needed to re- 
alize OSMPs, whereas beyond the critical bandwidth ratio of 
0.6 no OMST is possible within this model. 



mapped onto a self-consistent two-orbital Anderson im- 
purity modeliiSS with effective action: 



Jo So dTdT ' Em, dlAr)Q^ (r - r')d m<J {r') 



•f / ( fdr(n 1 +n 2 -2) 2 



(27) 



The hybridisations to the effective conduction bath are 
self-consistently related to the local interacting Green's 
functions G m through: 



g m (iu n ) 1 =iuj„- t 2 m G m (iu n ) 



(28) 



These equations are exact for an infinite-connectivity 
Bethe lattice (corresponding to semi-circular non- 
interacting d.o.s). Particle-hole symmetry with one elec- 
tron per site in each band has been assumed. We focus 
here on the paramagnetic solutions only. The DMFT 
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FIG. 11: (color online). Same phase diagram as in fig. [§]but 
in the U-J plane at T = and at /3ti = 40. 



equations will be solved in the following using both an 
exact diagonalisation (ED) and Quantum Monte-Carlo 
(QMC) technique. 



A. Exact diagonalisation study 

Within the adaptative exact-diagonalization 
methodiiSi, the effective conduction-electron bath 
is discretized using a finite number of orbitals N s . 
Hence, one considers the two-orbital Anderson impurity 
hamiltonian: 



N, 



Haim=^2 ^2 e i™a\ rn(J a lma + ^ ^ ^ m (d^ CT a ;TO(T + h.c) 

(29) 



ma l — l 



U 2 

+"2 (ni + h 2 - 2) 2 



The operators ai ma , a l ma describe the discretized 
conduction- bath degrees of freedom. The effec- 
tive parameters {e; m ,Vj m } have to be determined self- 
consistently, according to (|28|l . namely: 



JV, 

i=i 



\V h 



(30) 



The ED method becomes an asymptotically exact solver 
of the DMFT equations in the limit N s — > oo. In practice 
however, one can handle only a finite number of effective 
sites. For the case at hand, we used a T = Lanczos 
algorithm, with N s — 5 (i.e 5 effective sites per orbital). 
The self-consistency l|30[) is implemented on a Matsubara 
grid corresponding to a (fictitious) inverse temperature 
j3 (taken to be in practice in the range 200 to 500, which 
insures a good resolution on the low-energy physics). 
We monitor in particular the quasiparticle weights, ap- 
proximated as: Z m = [1 — lmY, m (iu!o) /luq] -1 (where 
w B =7r//3(2n + l)). 
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Our ED results for these quantities are displayed on 
Fig. El and compared to the slave-spin results, both as 
a function of U for fixed t2/ti and as a function of tijti 
for fixed U (inset). It is clear from this figure that, within 
the energy resolution which can be reached with N s = 5, 
an orbital-selective transition is indeed observed in the 
DMFT(ED) results when t%/ti is smaller than a critical 
value (close to 0.25), in remarkable agreement with the 
slave-spin mean-field. The quantitative value of the crit- 
ical coupling U/ti for the localisation of the wider band 
is overestimated, as usual in Gutzwiller-like schemes. 




FIG. 12: (color online). Quasiparticle residues at J = in 
DMFT(ED) for tijtx =0.1 (symbols). For comparison, the 
slave-spin mean-field approximation is also displayed (contin- 
uous line). Inset: same quantities as a function of t%jt\ at 
fixed U /t\ = 4.0. (The kinks in Z\ are associated with the 
criticality of Zi)- 

The phase diagram obtained with ED, as a function 
of ti/t 2 and U/t\ is displayed in Fig. El Good quali- 
tative agreement with the slave-spin mean field is found 
(Fig. El. We also studied the hysteresis properties by 
performing runs for increasing and decreasing values of 
U/t\. This results in two almost parallel transition lines 
on Fig. El on e corresponding to the disappearence of 
the metallic solution (obtained from a series of runs for 
increasing U) , the other corresponding to the loss of the 
insulating nature of the wide band (from a series of runs 
for decreasing U). In the region between these two lines, 
coexistence of two types of DMFT solutions is found: for 
small t-zjty, one of the solutions is orbital-selective Mott 
and the other is fully insulating, while for larger t<ijt\, 
one of the solutions is metallic and the other fully insu- 
lating. The actual thermodynamic transition is given by 
the crossing of the free energies of the two solutions. In 
contrast, no hysteresis has been found at the transition 
between the OSMP and metallic phases. This suggests 
that the transition from the metallic to the OSMP phase 
is of a very different nature than the Mott transition of 
the wider band. If only the gap closure is monitored, 
then only one transition is found at J = 0, in agree- 
ment with the symmetry argument of RefA As we shall 
demonstrate below, there is indeed strong evidence that 



the orbital-selective Mott phase at J — does not display 
a sharp gap in either orbitals. 
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FIG. 13: (color online). Phase diagram at J = and T = 
in ED-DMFT. The dashed line marks the disappearance 
of the Mott gap (downward runs). For U values around 5 
there is coexistence between an OSMP and - depending on 
the bandwidth ratio - an insulating or metallic phase. 



B. Low-energy nature of the orbital-selective Mott 
phase at J = 0, from ED and QMC 

In order to understand better the nature of the orbital- 
selective phase found at J = 0, we take a closer look at 
the local Green's function for each orbital, in each phase. 
In order to do this, we also solved the DMFT equations 
using the Quantum Monte-Carlo algorithm of Hirsch and 
FyeSS. The ED and QMC methods are quite comple- 
mentary. The former applies at T = but suffers from a 
limited energy- resolution due to the small value of N s , 
while the latter is limited to finite-temperature but can 
be made very precise by increasing the number of time 
slices and the number of Monte-Carlo sweeps (we used 
128 slices in imaginary time and up to 5 x 10 5 Monte 
Carlo sweeps in practice). Also, using QMC allows for a 
reconstruction of the spectral functions using a numeri- 
cal analytic continuation based on the maximum entropy 
algorithm. 

In Fig. El wc display the ED results for the local 
Green's functions on the Matsubara axis, for a small 
bandwidth ratio t^/ti and three different values of U 
corresponding to the insulating, metallic, and orbital- 
selective Mott phases. In the particle-hole symmet- 
ric case, the Green's functions are purely imaginary on 
the Matsubara axis and related to the spectral function 
A m {e) of each orbital by: 

ImG m («j) = -2w / + de A ^ e \ (31) 
Jo w + e 

When the spectral function A m (e) has a gap, the in- 
tegral in the right-hand side of i|31|) has no singularity 



11 



in the u> —* limit, and hence ImG m (iuj) a w at low 
frequency. Furthermore, ImG m (iw) has a minimum for 
lj of order A m /2, with A m the gap in the m-th or- 
bital (as can be seen by replacing A m (e) by the sim- 
plified form l/2[S(e - A m /2) + 6(e + A m /2)], yielding 
lmG m {iuj) ~ -2w/(w 2 + A^/4)). This is fully consistent 
with the ED results in the upper panel of Fig. 1141 corre- 
sponding to a large value of U/tx = 7, with both orbitals 
insulating and having a gap of order U. 

In the lower panel of Fig. ^] ED results are displayed 
for U jt\ — 2, when both orbitals are metallic. In this 
case, the low- frequency limit of (|^T|) yields: lmG m (iuj — > 
0) = — wA m (Q). This is consistent with the numerical 
results, which also show that the Luttinger theorem is 
obeyed for both bands: TrA m (0) — l/t m . 

The central panel of Fig. ITU displays the ED results for 
an intermediate coupling, corresponding to the orbital- 
selective phase. One sees that the wider band has metal- 
lic behaviour, with A\{uS) still reaching the Luttinger 
value at u = 0. In contrast, lmG2(iu>) appears to vanish 
as lj — > 0, within the energy resolution of ED. However, 
in striking contrast to the upper panel (insulating phase) , 
the minimum in ImG2(iw) is at a very low frequency scale 
which is obviously not given by U. This strongly suggests 
that A2(oj) displays low-energy peaks very close to lj = 0, 
and may even have spectral weight down to arbitrary low 
frequency. 

Figure ^| compares the ED and QMC results for 
ImGi j 2(«w) in the strongly anisotropic case t<ijt\ = 0.1 
for a relatively small Coulomb interaction Ujt\ = 1.6, 
easier to study with QMC. These parameters also cor- 
respond to the orbital-selective phase (Fig. I13|) . Very 
good agreement between the two methods is found, con- 
firming the above analysis (and confirming also the Lut- 
tinger value for the wider band with greater accuracy 
than in ED). The corresponding spectral functions ob- 
tained by the maximum-entropy method are displayed 
on Fig. The spectral function of the broader band is 
only slightly modified as compared to the non-interacting 
d.o.s. Small shoulders are visible, at the position of the 
lower and upper Hubbard bands, the Luttinger theorem 
is obeyed, and some of the spectral weight is transfered to 
higher energies as expected. The narrow-band however is 
obviously in a strong-coupling regime, with well-marked 
upper and lower Hubbard bands. The most striking fea- 
tures however are the two narrow peaks at low-frequency, 
which can be interpreted either as a split quasi-particle 
resonance or as the sign of a pseudo-gap (partially filled 
by thermal excitations since the QMC calculation is for 
T/tx = 1/40). 

Hence, the general conclusion of this analysis is that 
the orbital-selective phase found for J = at small 
enough ti jt\ is not a conventional Mott phase in which 
the ( "localized" ) orbital with narrower bandwidth would 
display a sharp gap. Instead, two narrow peaks exist 
near u> = and finite spectral weight is found down 
to low-energy. Our numerical data are consistent with 
a pseudo-gap behaviour, but a precise characterization 
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FIG. 14: Imaginary part of the Green functions of the two 
bands at low energy for <2/ti = 0.16. Upper panel: U/ti = 
7.0, both bands are insulating. Central panel: U jt\ = 4.0, 
Orbital Selective Mott phase. Lower panel: U/ti — 2.0, both 
bands metallic. The dashed curves correspond to the orbital 
with narrower bandwidth. All energy scales are in units of ti. 



of the low-energy nature of the phase will require fur- 
ther effort, using highly precise techniques at low-energy 
such as the numerical renormalization group. This is 
left for future work. Such a study should also clarify in 
which precise sense the narrower band is "localized" in 
this phase, and whether the orbital-selective transition is 
a true phase transition or rather a sharp crossover. 
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FIG. 15: (color online). Imaginary parts of the Green's 
functions in Matsubara space for fo/ti = 0.1, U/ti = 1.6. 
Solid and dashed lines represent the exact diagonalization 
results for the wide and narrow bands respectively. Circles 
and squares represent QMC data for the same quantities at 
/3ti = 40. 




FIG. 16: (color online). Spectral function for i2/ti = 0.1, 
U/ti = 1.6 at pti = 40 in QMC-DMFT. The solid and dashed 
red lines denote the narrow and wide bands respectively. The 
DOS of the non-interacting system (wide band only) is given 
for comparison. 



VI. INSTABILITY OF THE 
ORBITAL-SELECTIVE MOTT PHASE UPON 
INTER-BAND HYBRIDIZATION 

To what extent an orbital-selective Mott phase may 
occur in practice depends on its stability with respect to 
perturbations. This is an important issue in view of the 
fact that the hamiltonian considered in this paper has a 
rather high degree of symmetry. In this last section, we 
consider the effect of an hybridization between the two 
orbitals (also recently considered in Ref^), i.e of a local 
non-diagonal term: 



H, 



hyb 



di2<rdilcr) 



(32) 



We note that this term could be eliminated by diagonal- 
izing the non-interacting hamiltonian. However, in the 
new basis, the interaction terms will be modified: terms 
will be generated which will have the same physical ef- 
fect than a hybridisation (and will involve non-local con- 
tributions in general). Indeed, as emphasized in Ref. 23 , 
the existence of an OSMT is a basis- independent issue. 
In a general two-band model, a Mott transition is sig- 
naled, when approached from the metallic side, by a low- 
frequency singularity in loI—Y,(lo) = Z~ 1 lu+- ■ ■ , where S 
and Z are the self-energy and quasiparticle- weight matri- 
ces, respectively. An OSMT is characterized by Z having 
one zero-eigenvalue while the other one remains finite. 
Being associated with the rank of the Z-matrix, it is a 
basis-independent notion. Our choice of basis is such that 
the interaction terms have the form specified above. 



A. Physical considerations 

Some statements about the effect of a finite hybridiza- 
tion can be made on general physical grounds (focusing 
for simplicity on the J = case). First, for small values 
of U when both orbitals are itinerant, it is obvious that a 
hybridization will not change qualitatively the low-energy 
nature of the metallic phase. Also, at very large U, when 
both orbitals are localized and a gap exists in both or- 
bital sectors, we expect the presence of a gap to be a 
robust feature which persists in the presence of V . How- 
ever, it is also clear that introducing V into this gapped 
insulator allows the local moments formed in the Mott 
insulating state at V = to screen each other. This 
occurs through the formation of on-site bonding and an- 
tibonding "molecular" levels mixing the two orbitals. As 
a result, the local- moment Mott insulator is expected to 
be replaced by a Mott insulator in which local singlets are 
formed on each site. The intermediate U regime, in which 
the system is in the OSMP at V — 0, is more delicate. Be- 
cause orbital 1 is itinerant, a Kondo screening process can 
take place, which will screen the local moment (formed 
by orbital 2 when V = 0). The resulting state can a 
priori be either a heavy-fermion metallic state involving 
quasiparticles with a large effective mass, or (because we 
are considering the half-filled case), a Kondo insulating 
state in which a gap is formed in the low-energy quasi- 
particle spectrum. Below, we study this question using 
the slave-spin approach and find that both phases can 
be obtained, depending on the value of U and V. A 
larger V favours the opening of a gap, as expected. Us- 
ing a general low-frequency analysis, we also demonstrate 
that, for small values of V, the heavy-fermion metallic 
state, and not the Kondo insulating phase, is induced. 
We note, finally, that both the Kondo-insulating phase, 
and the phase obtained at large U in which on-site local 
moments in different orbitals screen each other, have a 
singlet ground-state. It is therefore not obvious a priori 
whether these two phases are continuously connected^, 
as found in2S*2i for a related (but different) model, or 
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whether a phase transition between them can exist. We 
return to this point at the end of this section. 

In any event, these physical arguments imply that the 
orbital-selective Mott phase is unstable with respect to 
the introduction of a non-zero hybridization, at zero- 
temperature. Of course, for temperatures above the 
Kondo scale, the physics of the OSMP can be recovered. 
This is an important point in view of the possible exper- 
imental relevance of the OSMP. 



B. Slave-spin mean- field study 

We now turn to a more quantitative study of the effect 
of a finite hybridization (|32(l , using the slave-spin mean- 
field approximation. Our aim is not to establish a full 
phase diagram for all values of the parameters U, J and 
V, but rather to investigate whether this approach does 
support the physical expectations discussed above. A 
more extensive investigation will be presented in a future 
publication. 

In Fig. El we display the quasiparticle residues Z\ (for 
the broad band) and Z 2 (for the narrow band), as a func- 
tion of U and for increasing values of the hybridization 
V. It is seen that, starting from the OSMP in which 
Z\ 7^ and Z-i = for V — 0, one obtains either a 
phase in which both Z\ and Z2 are non-zero (i.e Kondo 
screening takes place), or an insulating phase in which 
Z\ = Z-2 = 0. This demonstrates that the hybridization 
is indeed a singular perturbation on both the J = and 
finite- J orbital-selective Mott phases, in agreement with 
the physical arguments above. 

In order to check whether the Kondo-screened phase is 
a (heavy-fermion) metal or whether it is gapped (Kondo 
insulator), we have plotted in Fig ^| the band gap of 
the auxiliary quasiparticlcs found within slave spin mean- 
field theory (the plot is for J = 0, the J ^ case being 
similar). Note that due to the factorization of the Green's 
function 

(^4(0)) = 4(S x (0) mry S x (r) mry )(f mry (T)fl a (0)) 

(33) 

the physical gap in the insulating phase will be of the or- 
der U. For small and intermediate values of V, the phase 
with Z\ 7^ and Z 2 ^ is metallic (gapless) (see Fig ITS)) . 
The orbital-selective Mott phase is replaced by a heavy- 
fermion regime, as shown in Fig. 1171 in which the orbital 
with narrower bandwidth acquires a very large effective 
mass (corresponding to a very low quasiparticle coher- 
ence scale Z2). This is also in qualitative agreement with 
our recent study of the periodic Anderson model with 
direct f-electron hopping^. Only beyond a critical value 
of V is a gapped Kondo insulator found. As discussed 
below, this can in fact be proven generally, beyond the 
mean-field approximation used here. 

At small to intermediate values of V, the two Mott 
transitions associated with the OSMP at V = are there- 
fore replaced by a single non-selective transition from a 



(heavy-fermion) metal to an insulator as U is increased. 
Within the slave-spin mean-field, this unique metal to in- 
sulator transition is found to be first-order, and to occur 
at a critical value of U which lies in between the criti- 
cal interactions of the metal-OSMP and OSMP-insulator 
transitions. 

For larger values of V, the OSMP phase is replaced by 
a Kondo- insulating phase with Z\ , Z2 7^ but a finite 
quasiparticle band- gap (lower plot in Fig. IT7|) . Within 
slave-spin mean-field, a phase transition takes place as 
U is increased, towards another insulating phase with 
Z\ = Z2 = (corresponding to the fact that the Kondo 
effect does not take place when V is turned on starting 
from a Mott phase for both orbitals, with a large gap). 
As pointed out above however, this large- U insulating 
phase also has a singlet ground-state however, due to 
inter-orbital screening. Whether this phase transition is 
an artefact of the slave-spin mean-field or whether it is 
indeed present in a more accurate DMFT treatment is a 
question to which wc shall return below. 




12 3 4 

U/t, 




01234567 

u/t, 

FIG. 17: (color online). Quasiparticle residues Z\ and Z2 
within slave spins MFT, for finite V. Top: t 2 /ti = 0.5, J = 
0.2517 for (from right to left in the upper manifold -wide band 
-, from left to right in the lower manifold - narrow band): 
V = 0,(OSMT system),V/ti = 0.1,0.15,0.2,0.3. Bottom: 
tz/ti — 0.1, J = for (from right to left in the upper manifold 
-wide band -, from left to right in the lower manifold - narrow 
band) V = 0,(OSMT system), V/ti = 0.05,0.1,0.2. Insets 
show the same graph in log scale. 
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From this expression, it is easily seen that a gap is present 
whenever: 



0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 

U 

FIG. 18: (color online). Gap amplitude of the auxiliary 
fermions (see text) and quasiparticle residues within slave 
spins MFT, for two values of V. Top: t 2 /ti = 0.1, J = 
for V/t\ =0.1 (upper panel), and V/ti = 0.6 (lower panel). 
The corresponding critical V for the non-interacting system 
is \JDiDi ~ 0.63 in this case. 



C. General low- frequency analysis 

In order to understand better the nature of both the 
metallic and the insulating phases to which the OSMP is 
driven for V =/= 0, we perform here a low-energy analysis 
in terms of a renormalized quasiparticle band-structure 
(see also2£). We focus on the case in which the Kondo 
effect does take place as V is turned on, resulting in both 
quasiparticle residues being finite. We thus keep the fol- 
lowing terms in the low-frequency expansion of the self- 
energy: 



Eii(w) = w(l-1/Zi) + 
£ 22 H=lj(1-1/Z 2 ) + 
E 12 (w) = Eia (()) + ••• 



(34) 



The low-energy quasiparticle band structure is then given 
by: 

[u - Ziei (k)] [u - Z 2 e 2 (k)\ - Z x Z 2 [V + E 12 (0)] 2 = (35) 



(36) 



In this expression, V c g is the effective hybridisation 
between low-energy quasiparticles, and -Dieff, -D2eff the 
renormalised quasiparticle half-bandwidths, given by: 



V eS = v / ^i^ 2 [V-+E 12 (Q)] , D lcS = Z X D X , D 2cff = Z 2 D 2 

(37) 

Remarkably, Z\ and Z 2 drop out from 1)36(1 , and the con- 
dition for a quasiparticle bandgap therefore reads: 



V + E 12 (0) > y/DtDa 



(38) 



If this condition is satisfied (and Z\, Z 2 ^ 0), one has 
a Kondo insulator phase, with a quasiparticle band gap 
given by: 



(D leS - D 2eS y + W C 2 S - (D lcff + D 2cS ) (39) 



In the opposite heavy-fermion (gapless) metallic 

phase is formed. 

From l|38l) . it is seen that it is the off-diagonal com- 
ponent of the self-energy (induced by V) which plays 
the key role in deciding whether a gap opens or not. 
From the same equation, it is also clear that a band- 
gap cannot open for arbitrary small V . Indeed, at small 
V, E12 grows at most proportionally to V, and hence 
the left-hand side of l|38|) is small, while the right-hand 
side is finite (note that the r.h.s involves the bare band- 
widths). Hence, a critical value of V is required to open 
a band gap and enter the Kondo insulating phase, start- 
ing from the OSMP. Note that this analysis is general 
and relies only on Fermi-liquid considerations, indepen- 
dently of the specific method used to solve the model. 
Note also that the situation is very different when one of 
the bare bands is dispersionless (D 2 =0), as is the case 
of the periodic Anderson models considered inS^i. In 
this case, it is clear from (|38|) that an arbitrarily small 
V induces the Kondo insulating state. The formation 
of a heavy-fermion metallic state induced by a non-zero 
hybridization in a generalized periodic Anderson model 
with D 2 ytz (direct f-electron hopping) was also investi- 
gated in our recent work2i 

The slave-spin mean field results presented above can 
be placed in the context of this general low-frequency 
analysis. Within this approach, one has: 

Z 1 = 4(Sf) 2 , Z 2 = 4(5f) 2 , V cS = 4(5f Sf) V (40) 

It is important to realize that this approach does provide 
an off-diagonal component of the physical electron self- 
energy, which, using H37fl . is given by: 



E 12 (0) <Sf5f)-(5f)(5f) 



V 



(Sf){S$) 



(41) 
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Note that in more conventional slave-boson approaches, 
the physical electron operators are related to the quasi- 
particles by a slave-boson condensation amplitude which 
is a c-number {S m = V Z m f} n ). As a result, V c s — 
\fZ1Z2V , and no off-diagonal component of the self- 
energy is present. In such slave-boson approaches, the 
condition for the presence of a quasiparticle band gap 
is therefore entirely unrenormalized by interactions and 
reads: V f \J D\D% > 1. This is an oversimplification 
which is not present in the slave-spin approach. There, 
the criterion for the opening of a gap reads: 



V 



> 



(42) 



However, despite the renormalisation of the hybridisa- 
tion by the off-diagonal self-energy, we find that in prac- 
tice this criterion is very close to the non-interacting one, 
which gives a reasonable approximation of the critical hy- 
bridization necessary to open a band gap. We also note 
that, within the slave-spin approximation, the insulat- 
ing gap takes its bare value 2V (corresponding to the 
splitting between the bonding and antibonding molecu- 
lar levels) as soon as Z\ = Z 2 = 0, which is certainly an 
over-simplification. 

Finally, we comment on the phase transitions between 
the different phases induced by a non-zero hybridization. 
At small V, it is clear that there must exist a phase tran- 
sition between the metallic (gapless) heavy-fermion phase 
and the insulating (gapped) phase, as U is increased (as 
indeed seen on Fig. 1181 top panel) . The situation is less 
clear at larger V. There, we always expect an insulator 
with a singlet ground-state. However, the mechanism be- 
hind this singlet formation is rather different at smaller 
and larger values of U . In the former case, Kondo screen- 
ing dominates and the singlet is formed by screening the 
local moment in orbital 2 by the electrons in orbital 1. 
At larger U, the Kondo coupling is smaller, and it is 
the formation of an inter-orbital molecular bonding level 
(due to V) which is responsible for the singlet forma- 
tion. Within slave-spin mean theory, we find a phase 
transition between these two kinds of insulator as U is 
increased, the latter one being signalled by Z\ = Z2 = 
(Fig. [^1 bottom panel). However, it is not clear whether 
this phase transition is a real feature or an artefact of 
the slave-spin approximation. In—, Koga et al. recently 
suggested that these two phases should be adiabatically 
connected, as found also in the study of the periodic An- 
derson model with correlated conduction electron s 26 ! 27 
(note however that the present model is different, in that 
both bands have a dispersion and that an inter-orbital 
interaction is present). The problem is qualitatively sim- 
ilar (but not equivalent, because of the inter-orbital inter- 
action) to the two-impurity Kondo problem 2 ^, in which 
the inter-impurity singlet (RKKY) fixed points and the 
Kondo- singlet fixed points are in general adiabatically 
connected (except in the special case of particle-hole sym- 
metry, where a phase transition does occur) . Note that it 
is well knownSi that for this latter model, slave boson ap- 



proximations do lead to spurious first-order transitions. 
A full answer to this question is beyond the scope of this 
paper, and is left for a future investigation, together with 
a complete phase diagram as a function of interaction 
strength and hybridization. 



VII. CONCLUSION 

In this article, we have studied whether the Mott tran- 
sition of a half-filled, two-orbital Hubbard model with un- 
equal bandwidths occurs simultaneously for both bands 
or whether it is a two-stage process in which the orbital 
with narrower bandwidth localizes first (giving rise to 
an intermediate 'orbital-selective' Mott phase). In or- 
der to study this question, we have used two techniques. 
The first is a mean-field theory based on a new repre- 
sentation of fermion operators in terms of slave quantum 
spins. This method is similar in spirit to the Gutzwiller 
approximation, and the slave-spin representation has a 
rather wide range of applicability to multi-orbital mod- 
els. The second method is dynamical mean-field theory, 
using exact diagonalization and Quantum Monte-Carlo 
solvers. 

The results of the slave-spin mean-field confirms sev- 
eral aspects of previous studie s 9 ! 11 , and in particular the 
possibility of an orbital-selective Mott transition. How- 
ever, some of the conclusions differ from those of previous 
work. Specifically, the slave-spin approximation suggests 
that a critical value of the bandwidth ratio (t 2 jt\) c exists, 
such that the Mott transition is orbital-selective for ar- 
bitrary value of the Coulomb exchange (Hund coupling) 
J when t 2 /ti < (t 2 /ti) c . When t 2 /*i > [te/tijc, J has 
to be larger than a finite threshold for an OSMT to take 
place. This suggests that the existence of an OSMT is 
not simply related to the symmetry of the interaction 
term only. In particular, an intermediate phase is found 
for J = at small t 2 /ti. 

We have studied whether DMFT confirms these find- 
ings, and found that the main qualitative conclusions on 
the existence of the orbital-selective phase are indeed the 
same, but that the nature of the intermediate phase at 
J = is a rather subtle issue. Indeed, the narrow band 
does not have the properties of a gapped Mott insula- 
tor in this phase and displays finite spectral weight down 
to arbitrary low-energy. This is, for example, consistent 
with a pseudo-gap behaviour but requires further studies 
to be fully settled (using e.g low-energy techniques such 
as the numerical renormalization group) . 

We note also that our study emphasizes the key role 
of the exchange and (on-site) inter-orbital pair hopping 
terms in the Coulomb hamiltonian in stabilizing the 
orbital-selective phase, in agreement with Koga et alJi. 

Finally, we found that the orbital-selective Mott phase 
is generically unstable with respect to an inter-orbital hy- 
bridization V. In the presence of such a term, two possi- 
ble phases are obtained, depending on the strength of U 
and V. Either the narrow orbital acquires a large (but 
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finite) effective mass, corresponding to a heavy-fermion 
metallic state. Or the system is an insulator with a 
gap. This insulator differs from the Mott insulator at 
V = since it has a singlet ground-state. This is due to 
screening processes, involving both Kondo exchange and 
the formation of an on-site molecular (bonding) level. 
Whether one has in fact two different insulating phases 
separated by a phase transition (each phase being dom- 
inated by one of these screening processes) -as obtained 
by slave-spin mean-field-, or whether one has a simple 
crossover— is an open question which deserves further 
study. 

Of course, at intermediate temperature (above the 
quasiparticle coherence scale of the narrower band, but 
below that of the wider band), a physics similar to the 
orbital-selective Mott phase can be recovered even in the 
presence of a finite hybridization. This orbital-selective 
heavy-fermion state might be relevant to the physics of 
Ca2- :E Sr :I .Ru04. This is indeed supported by the re- 
cent angular magnetoresistance oscillations experiments 



of Balicas et al—. 
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